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Abstract 

We present a novel approach, based entirely on the gravitational potential, for studying the evolution 
of non-linear cosmological matter perturbations. Starting from the perturbed Einstein equations, 
we integrate out the non-relativistic degrees of freedom of the cosmic fluid and obtain a single 
closed equation for the gravitational potential. We then verify the validity of the new equation 
by comparing its approximate solutions to known results in the theory of non-linear cosmological 
perturbations. First, we show explicitly that the perturbative solution of our equation matches 
the standard perturbative solutions. Next, using the mean field approximation to the equation, 
we show that its solution reproduces in a simple way the exponential suppression of the non-linear 
propagator on small scales due to the velocity dispersion. Our approach can therefore reproduce the 
main features of the renormalized perturbation theory and (time)-renormalization group approaches 
to the study of non-linear cosmological perturbations. We conclude by a preliminary discussion of 
the nature of the full solutions of the equation and their significance. 
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1 Introduction 



The Largc-Scale Structure (LSS) of the universe grows from an initial nearly scale-invariant spectrum 
of Gaussian fluctuations due to a gravitational instability. The matter distribution of our universe 
today is well described on large scales by linear perturbation theory about a homogeneous and 
isotropic background. On length scales below about 10 Mpc the dynamics of matter is highly 
non-linear, so to describe its evolution in a qualitative way one has to resort numerical N-body 
simulations. 

At intermediate (quasi-linear) scales, the evolution of matter may be described analytically by 
extending standard perturbation theory (SPT) [1]. In SPT one defines a series solution to the fluid 
equations in powers of the initial density field. The n-th order term of this series grows as the n-th 
power of the scale factor a (for a pressureless fluid) which affects its convergence properties and limits 
its applicability. Qualitative comparisons to simulations have shown that the domain of applicability 
of second-order perturbation theory (in the linear power spectrum) is limited, at redshift z = 0, to 
wavenumbers of about 0.05 h/Mpc. 

The need for improving theoretical predictions for the next generation of very large galaxy 
surveys has spurred many efforts to go beyond SPT. Through future observations of matter at 
intermediate scales, such as the observation of the baryon acoustic oscillations at scales of the order 
of 10^ Mpc, one will be able to probe the expansion history of the Universe and of the nature of dark 
energy. The statistical weight of galaxy surveys is higher at smaller scales, increasing quadratically 
with decreasing scales. Since the aim is to get to percent accuracy in observations and theoretical 
calculations the need of extending the domain of validity of the theory of non-linear cosmological 
perturbations to smaller scales is clear. 

Renormalized perturbation theory (RPT) [2] reorganizes the perturbation expansion in terms of 
different fundamental objects, the so-called non-linear propagator and non-linear vertices, to improve 
the convergence of SPT. When the vertex is approximated by its lowest order (tree-level) form, then 
the matter power spectrum can be expressed as an expansion in terms of the non-linear propagator. 
The dominant contributions to the non-linear propagator are identified and summed explicitly in the 
large-momentum limit, resulting in an exponential suppression of the power at short distance scales 
due to the matter velocity dispersion. This is the significant feature of RPT: the emergence of an 
intrinsic ultraviolet cutoff which improves the convergence of the perturbation series at small distance 
scales. Then, by matching the behavior of the non-linear propagator at high momentum with the 
one-loop propagator which is the a good approximation to the true propagator at low momenta one 
obtains a non-perturbative prediction for the non-linear propagator. Substituting this propagator 
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in the first few diagrams of the reorganized expansion then gives a non-perturbative prediction for 
the power spectrum. The latter is the sum of two terms: the first term is exponentiaUy suppressed 
and proportional to the initial spectrum and the second term represents the power generated by 
mode-coupling at smaller scales whose time- and scale-dependence are dictated by the non-linear 
propagator. 

RPT relies on the standard equations for the density and velocity fields. An alternative approach 
is to use the fluid equations to derive evolution equations for the power spectrum and higher order 
correlators directly. For instance, the closure theory approach [3] approximates the three-point 
correlator by its leading order expression in SPT. The Renormalization Group (RG) method [4] and 
its variant, in which time is the flow parameter [5], offers an alternative, and to some extent, a more 
compact way of getting the same results of RPT. Truncating the RG equations at the level of the 
four-point correlator leads to a solution that corresponds to the summation of an inflnitc class of 
perturbative corrections and reproduces the RPT results. The RG method can be also be extended 
to include higher-order correlators [6], the effects of massive neutrinos [7] and modified to account 
for non-Gaussian initial conditions [8]. 

Additional methods to extend SPT have been proposed. The path-integral formulation of the 
Vlasov equation in terms of the distribution functions [9], in which statistical observables, such as 
the power spectrum, can be expressed as functional derivatives of a generating functional. The 
Lagrangian resummation theory (LRT) [10] reproduces the SPT power spectrum at the lowest non- 
trivial order and yields a non-perturbative prediction for the power spectrum that corresponds 
to resumming an infinite set of terms in SPT. The main difference between LRT and the RPT 
and RG methods is that in LRT the total power spectrum (including the mode-coupling piece) is 
exponentially suppressed at high momentum. 

The goal of this paper is to present a novel approach to study the non-linear evolution of matter 
perturbations which is entirely based on the gravitational potential. The motivation for our approach 
is the fact that the dimensionless gravitational potential is estimated to be smaller than about 10~^ 
on all cosmological scales of interest. The gravitational potential is of order unity only near black 
holes. For example, the gravitational potential in galaxies is a few times 10~^ and in clusters of 
galaxies it is of order a few times 10~^. Therefore, the gravitational potential is potentially a good 
order parameter for a perturbative expansion. It is certainly a better expansion parameter than 
the matter density contrast which is of order unity already in the quasi-linear regime. Of course, 
linear perturbation theory for the gravitational potential does not suffice since space derivatives of 
the gravitational potential can be large and comparable to the background energy density. Here the 
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fact that Einstein equations contain at most two space derivatives greatly restricts the number of 
terms that can appear in the full equation for the gravitational potential. 

The basic tool of our approach is a non-perturbative equation for the gravitational potential 
which we obtain by integrating out the degrees of freedom of the non-relativistic fluid. From the 
non-perturbative equation it is possible to reproduce in a few and simple steps the known SPT results. 
Furthermore, it is also possible to recover in a simple and physically transparent way the main results 
of the RPT and RG methods. In particular, the exponential suppression of the propagator in the 
high-momentum limit is recovered in our approach by solving our non-pcrturbative equation in the 
mean-field approximation. In this approximation the high-momentum modes of the perturbation 
"see" the velocity fluctuations as a modification to background and the nature of the suppression 
acquires a clear physical interpretation. 

Our equation provides a framework for evaluating the errors that one makes in various approx- 
imation schemes and allows to deduce some general conclusions about the evolution of non-linear 
cosmological perturbations even without a detailed analysis of its exact solutions. As we shall see, 
the non-perturbative equation indicates that the gravitational potential (and therefore the density 
contrast) can not grow indefinitely as time evolves. 

The paper is organized as follows. In section 2 we derive the non-perturbative equation for the 
gravitational potential. In Section 3 we show how the SPT results, up to third-order, are reproduced 
from our equation. In Section 4 we discuss the mean-field solution of the non-perturbative equation, 
reserving Section 5 to the comparison with the RPT and RG methods. In section 6 we discuss in 
a qualitative way the nature of the full solutions of our equation and their possible significance. 
Finally, Section 7 contains our conclusions and expectations for additional directions into which our 
analysis can be extended. 

2 The non-perturbative equation for the gravitational 
potential 

The goal of this section is to derive a non-perturbative equation for the gravitational potential in the 
presence of a fluid of non-relativistic dark matter particles. The extension to the more realistic case 
of the cosmological constant (or dark energy) dominance is straightforward and will be presented 
elsewhere. As we have mentioned in the introduction, to derive a closed equation for the gravitational 
potential we need to integrate out the fluid degrees of freedom, which we do using Einstein equations. 
Here we make use of the fact that the Einstein equations that the gravitational potential satisfies 
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contain at most two space derivatives, so if we wish to take them into account, it is enough to consider 
at most terms that are quadratic in the gravitational potential and contain two space derivatives. 
The rest of the terms are completely negligible on all relevant scales as explained in the introduction. 
Our starting point is the metric in the Poisson gauge 

ds^ = -a^(r) [-(1 + 2$) dr^ + lui drdo;' + ((1 - 2^')(5ij + xn) dxMx^] , (1) 

where a(j) is the scale factor as a function of the conformal time r, $ and ^ are the lapse function 
and the gravitational potential respectively, are the vector perturbations and Xij the tensor ones. 
As we are interested in subhorizon scales, we will retain first-order terms in * and $ and only the 
would-be second-order terms with two space derivatives. In particular, the difference — is 
second-order in so we need to keep V^(* — $). 

We do not assume that the vector perturbations vanish a;^ = 0, rather use the fact that the 
divergence of c<jj vanishes, V^Wj = 0. Similarly with the tensor modes which are transverse and 
traceless. This leads us to use the traceless-longitudinal projection operator 



= _ 3j_r_ (2) 



to determine $ in terms of * from the Einstein equations. This procedure should be compared to 
the lowest order procedure in which one simply uses the ^ 7^ j Einstein equation to determine that 
$ = We also use only the longitudinal part of the (iO) equations so the contribution of the vector 
perturbations drops. 

As we already mentioned, the energy-momentum tensor is that of a cold pressureless fluid of 
density p and three- velocity vector it, 

T^v = pu^,u^. (3) 

As usually done in the literature and to facilitate the comparisons, we work in the so-called single 
stream approximation and set to zero the stress tensor. The four- velocity vector u^j, = {uq,u) is 
normalized since in our units c = 1. Finally, we choose units for which SttG = = 1. With these 
assumptions, the relevant Einstein equations are (see, for instance, Ref. [11]) 



(00): m'^ + 2V'^m = pul, (4) 

(iO) : Vj*' + nVi^ = ^p uo Ui , (5) 

(ij) : (2^" + 6?^*' - 2 ^7^2 _ 2^L^ ^ + (v*)^ _ (* - $) j Sij 

-^)-2Vi'^Vj'i! = puiUj. (6) 
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For future reference we define the (matter) density contrast S = {p — ])) /p, where p = STi^ is the 
background density. The density contrast can be expressed in terms of the gravitational potential 
using the (00) equation (4), 



The traceless-longitudinal part of the (ij) equation is given by 



r'^ (ViV,(* - $) - 2Vi^Vj^) = V'^ (puiUj) . (8) 

We may evaluate explicitly some of the terms (and divide the whole equation by —2), 

V2(* _ $) + (V^r)2 _ 3__ (V.^rV,.*) = --pu^ + (pUiUj) . (9) 

We will use Eqs. (9) and (5) to express V^(^' — $) in terms of ^ only which will allow us to derive 
an equation for ^' that does not involve $. Despite appearances, Eq. (9) is completely local. In fact, 
all terms have at least two space derivatives acting on the fields. 
Let us now look at 1/6 of the trace of Eq. (6) 

+ sn-^' -(n^- 2—^ ^ + -(v^)2 - -v^ (^ - $) = -pu"^ . (10) 

\ a / 6 3 6 

The next step is to substitute Eq. (9) into Eq. (10) and obtain 

V' + sn^' - (n' - 2^) * + ^(v*)^ - ^ (v,w,*) = {puiuj) . (ii) 

Finally, we can use the (Oi) equation (5) to integrate out the three- velocity Ui from Eq. (11) 

^2 _ 2— j * + -(V*)2 - (ViWj*) 



= 2^ ( - ' ). (12) 

and the (00) equation (4) to integrate out the energy density p 

^"+3n^'-(n^ - 2^-] (v.^v,*) = / (v.^- + ^v.^)(v,^- + ^v,^) \ 

(13) 

This is our proposal for the full non-perturbative equation for the gravitational potential. We expect 
Eq. (13) to provide an accurate description of the non- linear evolution of cosmological perturbations 
on all relevant scales. Deviations form Eq. (13) will occur for situations in which our assumptions 
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about the matter content of the universe become inaccurate. In particular, we have assumed that 
matter is cold and so its intrinsic pressure is small, and that intrinsic dissipation is small. 

We will show shortly that Eq. (13) reproduces the standard perturbative results: that the gravita- 
tional potential is constant at the linear order and that it grows with time at second- and third-order 
(linearly and quadratic in the scale factor, respectively) . However, even without knowing the exact 
form of its full solution, Eq. (13) tells us that ^ cannot grow indefinitely. Indeed, when the den- 
sity contrast grows with time, that is when (2V^^'/3'H^) becomes increasingly large, the right-hand 
side of Eq. (13) gets more and more suppressed so ^ stops growing and most probably reaches an 
attractor solution corresponding to a density contrast that does not evolve in time any longer. 



3 Reconstructing perturbation theory 

We now verify that the perturbative solution of Eq. (13) reproduces the known perturbative solutions 
for the density contrast S and the velocity u. First, we linearize Eq. (13) and obtain the standard 
equation 

+ 2,%^' - 2 (n'^ - 2— j ^' = , (14) 
whose solution in matter domination does not evolve with time 

^{x,t) = -^i^{x). (15) 

For future use we recall that the resulting solution for u is 

2V^T r 

where in the second equality we have used the explicit form of 7i for matter domination, 71 = 2/r. 
The corresponding solution for 6 can be obtained from its relation to \I', Eq. (7), so 

-^L = ^r^V'^'L = ^a(r)V2*L • (17) 

Next we solve for ^ to second-order in perturbation theory. The idea is to substitute in the second- 
order terms the linear solution. Therefore, we take the non-perturbative Eq. (13), expand it to 
second-order in and take only terms that do not vanish on the constant leading order solution 
(15), 



-2—\ + -{v^f - ^^iVi^Vj^] =-^^(ViW,*j 



(18) 
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For matter domination it reduces to 



whose solution is 



5 V^V' / \ 1 , , 



To evaluate ^'2 in an explicit way it is easier to calculate V^*2) 



V2*2 = 3^r2 (^^V^V^^ (Vi*LV,*L) - ^V^ (V*l)')) . 



Thus, we need 



and 



Substituting into Eq. (21) we find 



14 \ 3 
In momentum space it gives 

(^2(^3, t) = -- 
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:A;|r^^'2(A;3,r) 



5 1 



kl + kl 2 (fei • fca) 



+ 



(27r)3 7 (27r)3 [7 ' 2^"" kfk^ 7 kfk^ 

X ^(fci + ^2 -fes) '^L(fci,r)(5L(fc2,r), 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



which reproduces the standard second-order kernel for the density contrast [1]. Using Eqs. (5) and 
(24), it is also straightforward to recover the kernel for the (divergence of the) velocity 6 = 'V ■ u. 
One expands Eq. (5) to second-order, using the time dependence of ^2 



= p(w2 + 5lWl). 



(26) 



Using the lowest order relations, the equation for the divergence of the velocity at second-order 
becomes 

(27) 



n 
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In momentum space, using the explicit solution (24) the final result is, 



02(^3, r) 


f d^ki 


f d^k2 


n 


J (2vr)3 J 


' (2vr)3 [ 



(28) 



X (5D(fel + fe2-fc3)5L(fel,'r)(5L(fe2,T), 



which agrees with the standard Newtonian kernel. Finally, we compare the third-order solution of 
our equation to the known solution of SPT. The third-order equation is given by 



14 V^V-'' / \ 4 V^V' 

^'^ + m^'^ = - ■ v^l) + — ( Vi^2V,*L I - 



3 VM'-^'-^y V- ^V.WlV^*lJ. (29) 

Exposing the explicit time dependence of the second-order solution *2 = t^*2 and solving for the 
r dependence, we find 



^ 36 



(V^'2- 



\ 14 V^V-'' { ~ \ 1 V^V^' / n \ " 



(30) 



In momentum space 

A r ^3 



*3(fc4) = ^/|;^ j ^h{kl,k2,ks)6D{k,-k,-k2~k3) *L(fcl)^L(fc2)^L(fc3). 

(31) 

The function /s is a sum of three terms 



(1) _ 



"Y^fci • (^2 + ^3) 



5 {k2 + ks) ■ k2ik2 + ks) ■ k3 1,, , / 
3 \k^Tl^^ 2^^-^^) 



(2) 



(3) 



I(fc4-fcl)(fc4-(fc2 + fc3)) 

3 

1 (A;4 • fci)(fc4 ■ ^2) kl 
~9 fc| ■ 



5(fc2 + fe3)-fc2(fc2+fe3)-fc3 1„ , , 

3 \k^Tk^^ 2^^^-^^^ 



(32) 



Using the relationship (17) between Si, and we obtain 
1 



S3ik4) 



^3(^1,^2,^3) 



-i:ky^3ik4,T) 
b 



/ / / ^:;^-^3(fcl,fc2,fc3)<5D(fc4-fcl-fe2-fc3)(^L(fcl)'^L(fe2)'^L(fe3), 

/3(A;i,fe2,fc3) 



"-1 "-2 "-3 

We have verified that this coincides with the third-order kernel for the density contrast [1]. We 
expect that the perturbative series of the solution of our equation is equal to the solution in SPT. 
This can be verified by comparing the recursion relation of the perturbation series resulting from 
our equation and the standard kernels of SPT [1]. We postpone the explicit verification of the all 
orders relationship between the two kernels to a future investigation. 



(33) 
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4 The mean-field approximation 



After verifying that the solution of non-perturbative equation (13) reproduces the perturbative 
results at first-, second-, and third-order, we turn to improve the solution in a significant way, using 
the mean-field approximation. This will also allow us to make contact with the current techniques 
of solving for the density contrast beyond standard perturbation theory. 

We are interested mainly in the dynamics at small length scales for which linear perturbation 
theory for the density contrast is not valid. The idea is to use this separation of scales to average 
over the slowly varying (in space) large scale velocity fields and produce an improved linear equation 
for the gravitational potential which will take into account the effective dissipation that these large 
scale velocity fields induce. 

The simplest way to apply the mean-field approximation is to linearize the equation by setting 
higher order terms to their lowest order classical solutions. Implementing this approximation is 
more complicated if we start directly with Eq. (13) because the perturbative orders are mixed there 
(see however below). The easiest way to correctly identify the mean-field equation is to start from 
Eq. (6). One neglects all the second-order terms on its left-hand side and sets the right-hand side 
to be equal to "pS times the average of w^. 



where the average is over the initial conditions ^'l which is assumed here to be a Gaussian random 
field. All other terms only contribute a small correction to the background equation. We may now 
perform the average and restrict to matter domination, setting tq = 1, 



Alternatively, one can start from Eq. (13) and separate the perturbative orders carefully. First, it 
is necessary to use the form of the second-order solution Eq. (20) to find that ~ 'HV^2- This 
is then substituted into Eq. (26) from which one finds that to second-order. 



When this expression is substituted back into Eq. (13) and reexpressed using Eq. (16) in terms of 
the linear order velocity wl, one obtains after averaging exactly the mean-field equation Eq. (35). 
As we have discussed above, the mean-field approximation relies on scale separation. To under- 




(34) 




(35) 




(36) 
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stand this let us express the non-hnear term in momentum space, 



(V*w)^(ro) = - J J J -^q,^g,Uq,-Uq,6D{k-qi-q2-q3) 

d^q2 f d^qs 2 ,t, /o'7^ 

(2^ J (2^ \k-l2- Q3\ *fe-q2-93 ■"92 • Wqs • (37) 

If we now restrict the integral to the region for which k ^ q2,q3 (the "large-momentum" limit) then 

(V^* u^), (ro) . -kH, / / ll ' 
which in the mean-field approximation gets replaced by 

-^'^"i / 1^ / 1^ ■ ^-^) = -^'^^ / (Cf^ ■ = -k'^>.{uHro)) . (39) 

The operation of averaging is justified by the fact that the high-frequency mode '^k "sees" the 
other slower modes as a approximate constant background which in turn is equivalent to a spatial 
averaging (and to an ensemble averaging due to the ergodic theorem). Notice that the same result 
can be obtained starting from Eq. (11). In such a case the term in the right-hand side of that 
equation is the only one retained and gives in momentum space and in the high-momentum limit 
(at T = To) 

4 / /" d^g'2 f d^qs , ^ ^ \ 4 /" Sq 



97^23 

= -^k''^„{u\To)) . (40) 

The quantity = ^{u'^^tq)) is related in a simple way to a quantity that is widely used in the 
literature, 

crl = \<^l (41) 

where 

a^^lJ^Phq), (42) 

with Pj" being the linear power spectrum of the density contrast. To facilitate the comparison of 
our results to the literature we will use in what follows. 

Substituting the explicit time-dependence of H and going to momentum space we find from 
Eq. (35) 

K + Ik = -^k'^ay^k ■ (43) 
11 



To solve Eq. (43) we need to impose boundary conditions. We choose to impose boundary conditions 
that assume that the solution is exactly equal to the constant perturbative solution until time tq, 
so ^'fc(To) = ^'l and \I'^(to) = 0. Then for times r > tq the non-linear effects are "turned-on" and 
start to affect the solution. The solution of Eq. (43) with the specified boundary conditions is the 
confluent hypergeometric function, 

'9 1 



(44) 



*fe(r) = ^'fe(ro) o-'^i 
and can be expressed in terms of a Bessel function 

N-5/8 /c 

r^V2(r4 -4)] t(±\ [2fca,(r2 - r^)] . (45) 



The solution starts at unity for r = tq, decays exponentially for a while and then oscillates within 
an envelope that decays as r~^/^ with a periodicity that increases as r^. The relevant information 
about the solution is essentially contained in the solution of the first-order equation 

m^a' = 4V^* alr^ , (46) 

which is 

(*MF)fc (r) = ^fe(ro) exp ^-h^alir' - Tq^)) (47) 

in the region in which the argument of the exponential is less than about unity. The reason why the 
solution (47) suffices is that the validity of the mean-field solution itself is restricted to times r and 
wavenumbers k such that the exponent is less than or order unity. This can be seen by looking at 
Eq. (5). To derive the mean-field approximation we have used the lowest order solution for which 
= 0. However, to mean-field solution itself is time dependent, so our estimate of the average {v?) 
will no longer be accurate when the solution deviates substantially from a constant. Let us compare 
^' to %^ for the mean-field solution (47) when r » tq, 

n-^ = (48) 

r 

They become comparable when ^k'^a^r'^ ^ 1. This means that the exponential solution is completely 
sufficient in the region of validity of the solution. For instance, at redshift z = 0, the exponential 
decay is valid at least up to scales k ~ 0.25 h Mpc"^ and then it gets replaced by a power-law decay. 
A similar estimate for the region of validity of the solution can be obtained by comparing to 
respect to SH'S' in Eq. (43). 
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Even though we cannot trust the exponential decay or the solution at times and/or wavenumbers 
for which fe^c^r^ S> 1, we do expect the softening of the decay to persist for the true solution of 
the full non-perturbative equation. The reason is that we expect the velocities to decay as a result 
of the decreased gravitational potential, and when the velocities decay, the effective dissipation that 
they induce will be weaker, so the rate of decay of the gravitational potential will be smaller. 

The mean-field approximation is also not expected to be valid for combination of times and 
wavenumbers for which fc^c^r^ <^ 1. The reason is that the scale separation that we have invoked 
previously to perform the averaging operation is not valid. Under such conditions the expected 
change to the solution is that, instead of being described by an expansion of the form e~'2^^'^^'^* = 
1 — ^k'^a'^T^ + • • • , the solution will be modified to 1 — cik'^a^r^ + • • • where ci 7^ ^. The modified 
coefficient ci can be determined by a perturbative (one-loop) computation. 

It is possible, if one so wishes, to iteratively improve the mean-field approximation by using 
the solution at each given step to recalculate {u'^) for the next step, generating an expansion in 
exp (— /c^fT.yT'^/G) . Another improvement would be in matching the mean-field solution in a better 
way to the perturbative solution or combining both improvements. We will not pursue these possible 
improvements any further, leaving them to future research. 

4.1 The mean-field solution for the density contrast 

The density contrast and the gravitational potential are related via Eq. (4). One would naively think 
that to obtain the mean-field solution for S it is enough to substitute the exponential solution (47) 
into the relation between ^ and S and obtain the desired result. However, this is not true, as we 
now show. 

Erom Eq. (7) we find that the mean-field equation that S obeys is somewhat different the its 
counterpart Eq. (35), 

6" + n6' - = ^{u^)V^d . (49) 

In momentum space, substituting the explicit time dependence of matter domination, 

S'l + -S'k - ^^k = -^k^crySk . (50) 

The diflFerent coefficient of H rather than 37{ in front of the friction term in Eq. (49) causes the time 
dependence of the mean-field solution for S to change. The full solution of Eq. (50) can again be 
expressed in terms of Bessel functions similar to the solution (45). The full solution is subject to 
the same restrictions as before so the relevant information is contained in the approximate solution 
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that is obtained by multiplying the linear solution 5l which solves the equation 

<^L + ^'^L-:^'^L = o, (51) 

by the solution of the first-order equation 

-4' = -Ak^aySk , (52) 
r 

whose solution is 

4 = 4(To)e-^^''^'(^'-"o) . (53) 
So the final result for the mean-field solution is 

{SMF)k = ('5L)fee-|'=''^?(-'--o) (54) 

An alternative way to find the mean-field solution for 6 is by starting from the standard fluid 
equations for 5 and 6 = 'V ■ u, 

6' + {1 + S) 9 + u^'VaS = , (55) 

e' + 'He + ^n^s + v^u'^Vauf^ + u'^VaO = o. (56) 

We then take the time derivative of Eq. (55) 

S" + (1 -h 6)6' + 96' + n"Va5' + iu''yVa6 = , (57) 

and substitute Eq. (56) 

6" -{1 + 6) (no + ^n^6 + Vfiu'^VaU^ + u'^VaO^ + 66' + (u'^yVaS + u'^Va6' = . (58) 

Now we substitute again Eq. (55) in the form Tid = —'H6' — T-L66 — T-Lu"'Va6 and find after some 
straightforward algebraic manipulations, 

6" + 'H6'-^nH 

-6u''Va6 + {Hu" + (n")') Va6 -{6 + n"Va) [(1 + 6)6 + u"Va(^] = . (59) 

To identify the mean-field solution we now go to the "high-momentum" limit, which in this case 
means taking the highest power of gradients acting on 6; the equation simplifies enormously. The 
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leading term comes from the second derivative term u^V aU^^ and one takes only the terms 
where the gradients act on J, 

5" + ns' - ^n^d - u'^u^VaVpS = o . (eo) 

Solving this equation in the mean-field approximation, 

5" + nd' - ^n^S = {u''u^)VaVp5 = \ {u^} . (61) 

which is the same equation as we have found for S previously, Eq. (49). The solutions of the two 
identical equations are, of course, the same. 

5 Comparison to the RPT and RG methods 

The exponential suppression that we have obtained in the mean-field approximation matches per- 
fectly the exponential suppression obtained for the non-linear propagator (to be discussed shortly) 
within the RPT proposed in Ref. [2] and the RG method put forward in Refs. [4, 5]. We would 
like now to show that this is not a coincidence, but that the mean-field approximation applied to 
our non-perturbative Eq. (13) captures the essential physics in the large-momentum limit described 
within the RPT and RG methods. 

The two methods introduce an object that is called the "non-linear propagator" which replaces 
the standard perturbative propagator (the Green's function) and encodes the information about 
the initial conditions, thus significantly improving over the standard perturbation theory. In the 
diagrammatic approach of RPT the non-linear propagator is obtained through the resummation of 
an infinite class of diagrams, while in the RG method it is evaluated by solving renormalization 
group equations obtained by introducing a running ultraviolet cut-off. 

Let us recall the definition of the non-linear propagator, 

G^ik.r)=,Mr)^t(fi^). (62) 

where the indices a, 6 = 1,2 correspond to the different type of perturbations (pi = 6 and 02 = 0, 
gab is the linear propagator, denotes the n'th order expansion of the full solution ^ in powers 
of the linear perturbation 0i^2 and the average is an ensemble average over the initial conditions. 

In the explicit calculations the non-linear propagator is calculated to third-order ( "one- loop" ) in 
the small momentum limit, while in the large- momentum limit it was found to have an exponential 
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form after resumming an infinite set of diagrams. In particular, if 

*^"^= / -^2i...a„(felr-- ,kn,T)SD{k-ki M -^a, (^l) • • • <^a„ (^n) (63) 

is the series expansion in terms of the hnear 0's and of the standard perturbation theory kernels 
J^aal—anJ the non-linear propagator may be written as the infinite sum (corresponding to an infinite 
number of diagrams) [2] 

oo „ 

Ga6(fc, T)=5a6(r) + ^(2n + l)!! / J^S,''^'\k,k,,-ki,- ■ ■ ,kn,-kn,T)Pl{ki) ■ ■ ■ Pl{kn) , (64) 
n=l 

with 

"^ab {k,ki, —ki, • • • , kji, —kji, r) = -^afeai---a2n (^' ^1' ' ' ' ' km ~kn, T)Ua-^ ■ ■ ■ Ua^^ , = (1,1). 

(65) 

These expressions considerably simplify in the large-momentum limit for which 

/<5$i^"+^^(fc,T)\ , ^(a(r) -1)2" r 3 . fe2n /r/ x r / ^ ■ Qi 

\ XA (1 \ / =9ab{r) -— / d g2n---d qi X2n ■ ■ ■ Xi {di,{q2n) ■ ■ ■ OL{qi)) , Xi = , 

\ d(j)b{k) I {2n)\ J q2n---qi kqi 

(66) 

such that 

GabiKr) = gab{r) + gab{r)f2 [-k'crMr) - 1)']'" 



n=l 



= 5„(.(T)exp {-k'aiia - 1)^2) , large fe-limit. (67) 

This exponential form of the non-linear propagator in the large momentum limit can be matched 
to the full one-loop expression and then used to calculate the matter power spectrum (and other 
cosmological observables) for the whole range of momentum (large and small), 

Ps{k,T) = Gj{k,r)Pl'ik)+PMc{k,T), (68) 

where Gg = Gu + G12 and all contributions which are proportional to the initial spectrum of 
fluctuations are now included in the first term, which accounts for the remainder of the primordial 
power at a given scale after non-linear evolution, and thus has direct information on the linear power 
spectrum; by contrast, the second term Pmc represents the power generated by mode-coupling at 
smaller scales, and depends on the linear power at different scales than k through complicated 
convolutions. The important aspect of Pmc is that its time dependence is also dictated by the 
non-linear propagator. As a result, the convergence of the perturbation series for Pmc is drastically 
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improved. The RG technique of Ref. [4] offers an alternative, and to some extent, more compact 
way of getting the same result. 

Comparing the definition of the non-linear propagator to our discussion of the mean-field ap- 
proximation in momentum-space, Eqs .(37-39) for the gravitational potential and Eq. (54) for the 
density contrast, it is clear the in the large- momentum limit, the non-linear propagators are exactly 
equal to the mean-field solution of our equation. For instance, in our case the non-linear propagator 
for is 

Consider now the non-linear term in the equation in the leading order approximation when all fields 
are replaced by the lowest order fields ^'l, 

= y/l^/ll/ll^' ^^^^ (fe - 91 - 92 - 93) , (70) 

where we have substituted the expression for the first-order velocity ml from Eq. (16). In the large- 
momentum limit, taking the functional derivative with respect to ^'l(qi) and averaging over the 
remaining two variables, leads to the mean-field result for the gravitational potential. Of course, a 
similar procedure shows that the mean-field solution is equal to the non-linear propagator for 6 and 
9 in the large-momentum limit. 

When we solve Eq. (13) in the mean-field approximation we first (space) average some of the 
terms of the equation and then find the solution of the resulting linear equation. In the relevant 
region of parameter space the equation is first order in the time derivative and so the solution is an 
exponential of the averaged term. In the RPT method the order is different, as is made particulary 
clear when using the so-called a method [12]. There one solves a linear equation which is first 
order in the time derivative for the propagator (or other correlation functions). The solution of this 
equation is an exponential, which is then averaged over the initial conditions. The two procedures 
do not necessarily commute and so we do expect differences at some level. 

To calculate perturbative corrections with the non-linear propagator one needs to know its full 
momentum dependence and not just its large-momentum form. The proposal of the RPT and RG 
methods is to use the one-loop perturbative propagator in the low k limit, the exponential form in 
the large k limit and in intermediate scales use a specific interpolating function between the two 
limits. The use of the interpolating function is well motivated physically and matches the data and 
simulations well. Clearly, if we use the same approximations as proposed in RPT and RG, then 
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the results that we obtain for the cosmological observables will be identical. The final form (68) for 
the matter power-spectrum in our approach will correspond to the sum of the homogeneous non- 
perturbative equation in the mean-field approximation (corresponding to the first term in Eq. (68)) 
and the particular solution induced by the other non-linear terms of the equation (corresponding 
to Pmc) Ei'iid containing the mean- field exponentially-suppressed propagator. This has the effect 
of taming the bad ultraviolet behaviour of the solution for ^ at, for instance, one- loop. This step 
exactly follows the proposals of the RPT and RG methods. 

The interpolation procedure that connects the high-momentum limit to the perturbative, low- 
momentum limit, is not based on first principles as arc the two limits. One advantage of our 
non-perturbative Eq. (13) is that it allows to estimate the deviations of the non-linear propagator 
from the exact Green's function. For example, we have seen that the exponential form of the solution 
is restricted to times and scales such that the argument of the exponent is not much larger than 
unity. This is related to the fact that the true solution is a solution of a second-order equation while 
the exponential solution solves a first-order equation. In the RPT or RG methods, the exponential 
form of non-linear propagator is derived in the high-momentum approximation. In this limit the 
equation that the non-linear propagator obeys becomes a first order equation. It would be interesting 
to understand in detail the connection between our approach and the RPT and RG approaches in 
this context. 

6 Qualitative analysis of non-perturbative solutions 

The most significant qualitative feature that we can attribute to the exact solutions of our equation 
is that under generic conditions they do not allow the density contrast to grow without limit. This 
leads us to expect that in the full solutions the density contrast will saturate asymptotically to 
an attractor solution which is independent on time and so can only be space-dependent. Since we 
are always using comoving coordinates, this expectation does not contradict the theoretically well 
established (and observed) fact that local densities as observed by static observers can grow and 
become very large. 

To understand this point in a more quantitative way, let us look for a self-similar solution 
[13] of Eq. (13). The advantage of the self-similar ansatz is that it reduces the partial differential 
equation (13) to an ordinary (non-linear) differential equation which is much easier to solve. Here 
we restrict the discussion to a matter dominated universe. Since there is no preferred scale in the 
dynamics of a self-gravitating pressureless perfect fluid, cosmological fields should scale with a self- 



18 



similarity variable given some appropriate initial conditions. If we express the gravitational potential 
'^{x,t) as 

^{x,T)=T^^I;{y), y = x/T\ (71) 

then ~ TT-^ ~ rT"2 and V* ~ t'''. To make each of the terms in Eq. (13) scale with the 
same power of time, the power-law coefficients 7 and v must obey the relation 

7 = 2(z.-l). (72) 

In this case the density contrast S = 2V^*/3^^ (See Eq. (7)) scales as a constant in time. The 
scaling of the velocity can be read off of Eq. (5), u ~ t""^ in complete agreement with the known 
scaling solution. 

To fix the form of the scaling solution completely, that is to determine the parameter u, addi- 
tional input is needed. For example, it can be fixed by matching the matter power spectrum resulting 
from the scaling solution to the linear one. However, this procedure is not as straightforward as 
it seems because the initial conditions typically do not obey the scaling relations that the scaling 
solution obeys. The gravitational potential is a random variable on which one may impose initial 
conditions only statistically, e.g. by specifying the linear power spectrum for Gaussian initial condi- 
tions. A particular realization of the ensemble of initial conditions will not obey spatial scaling [14]. 
Furthermore, the introduction of (time-independent) smoothing scales and/or cut-off scales breaks 
self-similarity. Nevertheless, studying self-similar solutions is useful to confirm our expectation that 
on non-linear scales the density contrast saturates to a limiting value. 

An additional solution to our equation can be found using a simple ansatz 

^{x,t) = ^o{xo) + ^(^0)^^-7^ . (73) 

The velocity u vanishes for this solution; indeed, the (iO) Einstein equation can be expressed as 

^(a(T)V^)' = pti, (74) 

so the chosen ansatz corresponds to a vanishing velocity. Since Eq. (13) is non-linear the parameter 
A in the ansatz (73) is determined by the equation itself. It is possible to check that if ^ = — 1 the 
solution satisfies exactly Eq. (13). When we calculate the value of 6 using Eq. (7) we find that it 
corresponds to a constant density contrast S = —1, which is nothing but a void. This, in some sense, 
is a trivial solution whose existence should have been anticipated. However, the way it comes out of 
the equation is quite non-trivial and reassures us about the correctness of Eq. (13) and strengthens 
our expectations that the density contrast cannot grow in time without bounds. 
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Assuming various other forms for the solution, one can find additional candidate attractor solu- 
tions. For example, another candidate solution can be found by using a wave-like ansatz 

^ = Aln(^B -x- . (75) 

Here B and tq are constants. To understand the ansatz (75) better, let us look at its one dimensional 
version (taking the argument of the logarithm as dimensionless) 

= AlrL{x - vt) . (76) 

Using Eq. (7) we find 

b [X — VT)'^ 

When (x — vt) becomes small the magnitude of the density contrast grows and in this case we can 
approximate 6 ^ p. Using the same approximation we find 

12v , , 

u = ^ , 78 

SO, for large r, u decreases towards zero. Combining Eqs. (77) and (78) we find that the right-hand 
side of Eq. (13) is subleading compared to others for large times. 

Further, scales as Xjix — vt)'^ while ^f' ~ ll{x — vt) and is therefore subleading. Collecting 
the relevant terms, Eq. (13) reduces to 

*"-i(5,*)2 = 0, (79) 

which fixes the value of A in Eq. (76), 

A = -2v'^ . (80) 
This gives the final result for the density contrast (and so for the density) 



3 (x — vtY" 



P--o7Z—:iy2- (81) 



The solution (81) is reminiscent of the form of the density contrast in the Zel'dovich approximation 
[15], 

1 1 1 

~ l-a(r)Ai(q) l-a{T)\2{q) I - a{T)Uq) ' ^^^^ 
where x is the Eulerian coordinate and q is the Lagrangian one. The sign of the eigenvalues Aj deter- 
mines the nature of the solution. Comparing the Zel'dovich approximation to the solution (81), we 
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see that this solution corresponds to two equal positive eigenvalues and a vanishing third eigenvalue, 
corresponding to a collapse to a one dimensional filament. 

The full significance of the solutions that we have just discussed and of other non-perturbative 
solutions to our equation needs, obviously, to be understood in much more detail. We expect to 
carry this line of investigations to completion in the future. 

7 Conclusions and perspectives 

In this paper we have derived a non-perturbative equation for the gravitational potential in the 
case that the cosmic fluid consists of non-relativistic pressureless particles. Our equation reproduces 
the known perturbative results (verified explicitly up to third order in perturbation theory) and 
provides an alternative derivation of the exponential suppression of the non-linear propagator in the 
high-momentum limit obtained recently with the RPT and RG methods. 

We expect that our current findings can be extended in various directions. First, a quantitative 
comparison with observations. In this paper we were not pursuing a full quantitative computation 
of cosmological observables. This quantitative analysis and its comparison with data is of utmost 
importance as it is the main motivation for the study of non-linear cosmological perturbations. We 
also expect to be able to extend our results to the cases of cosmological constant (or dark energy) 
dominated universe and non-Gaussian initial conditions. 

Wc expect to be able to itcrativcly improve the mean-field approximation by using the solution 
at each given step to recalculate {u^) for the next step, generating an expansion in cxp (— /c^cr^r^/O) . 
This will allow us to extend the region of validity of the solution to larger values of fc^cr^r^. Another 
improvement would be achieved by improving the matching of the mean-field solution to the pertur- 
bative solution which is valid for small values of /c^(T,yT^ or combining both type of improvements. 

Finding exact non-pcrturbativc solutions or approximate attractor solutions can help in improv- 
ing the solution towards the fully non-linear regime. We do not expect our equation to be valid in 
the truly non-linear regime where intrinsic dissipation effects lead to virialized structures. We do 
expect to be able to understand its limitations in a quantitative way and consequently the resulting 
limits on the precision of its solution in the quasi-linear regime. 
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